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摘要 :河北 大 海 陀 自 然 保 护 区 地 处 我 国 暧 温带 落叶 阔 叶 林 区 , 草 甸 植被 多 处 于 山地 顶部 ,对 外 界 干扰 敏感 。 为 了 人 研究 该 地 区 的 
草 旬 植被 变化 以 及 其 对 气候 变化 的 响应 ,收集 了 该 地 区 近 30 年 的 TM 遥感 影像 资料 ,1980 一 2015 年 的 气候 数据 ,包括 年 平均 气 
温 、7 月 份 平均 气温 、1 月 份 平 均 气 温 .年 降水 量 .6 一 8 月 平均 降水 量 等 多 个 气候 指标 ;以 及 大 海 陀 自然 保护 区 及 其 周边 地 区 的 
地 形 数据 等 ,用 滑动 平均 法 、M-K 检验 法 .相关 分 析 法 及 偏 相关 分 析 法 年 方法 分 析 了 大 海 陀 自 然 保护 区 草 馈 植被 NDVI 及 其 与 
气候 响应 。 结 果 表 明 :(1) 近 30 年 来 ,大 海 陀 自然 保护 区 草 旬 植被 NDVT 昌 先 上 升 后 下 降 的 趋势 ,在 2004 年 NDVI 达到 最 高 , 随 
后 逐渐 下 降 。(2) 大 海 辽 自然 保护 区 草 多 区 域 的 年 降水 量变 化 整体 表现 为 波动 循环 的 趋势 ,总 体 咯 有 上 升 ,但 没有 达到 显著 水 
平 ;年 均 温 变化 表现 为 上 升 趋势 , 且 达 到 了 极 显著 水 平 ,该 地 区 的 气温 十 升 趋势 主要 由 以 1 月 为 主 的 冬季 温度 升 高 而 引起 。 
(3) 草 甸 植 被 NDVI 与 年 均 温 成 显著 的 负 相 关 关 系 ,与 年 降水 量 的 关系 处 明显 。 

关键 词 : 大 海 陀 自 然 保护 区 ;气候 变化 ; 草 甸 植被 ; 归 一 化 差 值 植 被 指数 
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Abstract: Hebe: Dahaituo Nature Reserve is located in a warm, temperate, deciduous broad-leaved forest area in China. 
Meadow. vegetation , which are sensitive to outside interference, are mainly distributed on top of Dahaituo Mountain. To 
determine the change and response of vegetation to climate change in the Dahaituo Nature Reserve and surrounding area, we 
collected a series of data, including Thematic Mapper( TM ) remote sensing image data for the last 30 years and climate data 
from .1980 to 2015, including annual average temperature, average temperature in July, average temperature in January, 
annual-precipitation, average precipitation in June to August, and terrain data. We analysed the relationship between the 
change in the meadow vegetation Normalized Difference Vegetation Index ( NDVI) and climate change using the moving 


average method, the M-K method, correlational analyses, and partial correlation analyses. The results showed that: (1) 
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NDVI of the meadow vegetation in Dahaituo Nature Reserve showed a downward trend in the past 30 years, after rising, 
achieved the best growth conditions in 2003, then gradually decreased. (2) The change in annual precipitation varied, but 
generally increased slightly. Annual average temperature showed a significant upward trend, with the temperature rise 
mainly caused by increased winter (January) temperatures. (3) NDVI of the meadow vegetation and average annual 
temperature were negatively correlated, whereas NDVI of the meadow vegetation and annual precipitation tended to be 


positively associated. 


Key Words: Dahaituo Nature Reserve; climate change; meadow vegetation; NDVI 


全 球 气候 变 上 暧 引起 了 全 球 平均 气温 上 升降 雨 格局 的 改变 ,这 些 环 境 因 子 的 变化 又 以 不 同 程度 不 同方 式 
影响 着 地 上 植物 群落 结构 ,物候 变化 以 及 地 表 和 土壤 微 环境 "。 植 被 是 陆地 生态 系统 的 主体 之 二 ,其 生长 状 
况 与 其 所 处 环境 密切 相关 ,如 气候 水质 和 地 形 等 环境 因子 。 气 候 变化 改变 了 植物 的 生长 环境 ,从 而 影响 到 了 
植被 的 生长 状况 “”。 国 内 外 学 者 已 开展 了 大 量 关于 植被 动态 变化 的 研究 ,其 研究 方法 主要 为 以 下 两 种 :一 、 
传统 野外 调查 与 观测 ;二 遥感 植被 监测 。 传 统 野 外 调查 方法 有 利于 获取 多 种 植被 类 型 的 详细 生长 数据 ,其 研 
究 内 容 较为 丰富 ,但 调查 所 需 时 间 长 且 费 力 ,造成 经 济 成 本 较 高 , 且 地 理 范 围 和 时 间 序 列 有 限 ,不 同调 查 人 员 
收集 的 数据 类 型 和 可 靠 性 存在 差异 ,因此 传统 野外 调查 方法 很 难 应 用 于 监测 大 区 域 尺 度 上 的 植被 动态 变 
化 ”。 自 20 世纪 70 年 代 , 卫 星 遥 感 技术 的 迅速 发 展 ,为 解决 上 述 问题 提供 了 良好 的 契机 “ 。 琐 感 影像 具有 
覆盖 面积 大 ,实时 性 、 周 期 性 准确 可 靠 性 等 优点 ,能够 从 景观 太 度 上 来 定量 监测 十 地 利用 的 变化 ,为 解决 当前 
资源 与 环境 之 间 的 问题 提供 了 有 力 的 技术 支撑 ,已 被 广泛 应 用 于 大 区 域 尺度 的 植被 动态 变化 研究 中 ” 。 国 
内 外 学 者 对 植被 与 气候 变化 之 间 的 关系 研究 较 多 ,时间 尺度 为 目前 国际 篆 用 的 时 间 序 列 ,主要 基于 归 一 化 卷 
值 植 被 指数 (NDVI) 数 据 研究 不 同 植被 类 型 对 气候 变化 的 响应 “” 。 

草 甸 是 以 多 年 生 .中 生 草 本 植物 为 主体 的 群落 类 型 ”。 在 草 甸 植被 的 研究 中 ,国内 学 者 通常 以 NDVI 为 
表征 ,探讨 草 旬 植 被 的 变化 趋势 及 其 对 气候 变化 的 响应 机 制 , 以 掌握 不 同 地 区 草 旬 植 被 的 演变 趋势 “'， ,研究 
区 域 主 要 在 我 国 西部 典型 草原 草 旬 分 布 区 s 河北 大 海 陀 自然 保护 区 山地 草包 与 我 国 西 部 草 匈 地 区 相 比 每 年 
6 一 8 月 平均 降水 量 和 平均 气温 相对 较 高 5 在 这 种 水 分 和 温度 相对 充足 的 情况 下 , 草 旬 植 被 的 变化 及 其 对 气候 
变化 的 啊 应 如 何 尚 没有 太 多 研究 与 报道 。 本 研究 以 草 旬 植被 的 NDVI 为 特征 指标 ,人 研究 河北 大 海 陀 目 然 保 护 
区 山地 草 馈 植被 变化 以 及 其 对 气候 变化 的 响应 ,为 进一步 探讨 我 国 东部 地 区 山地 章 馈 植被 变化 其 对 气候 变化 
的 啊 应 机 制 提供 支撑 。 


1 研究 区 域 


河北 大 海 陀 国家 级 自然 保护 区 地 处 河北 省 赤城 县 西南 部 ,其 地 理 坐 标 为 115°42'57” 一 115°57’0”E, 40°32' 
14” 一 40°41'40"N ,总 面积 12634hm”"” ,主峰 为 海 陀 山 ,海拔 2241m ” ,是 河北 赤城 县 与 北京 市 的 分 界 岭 “ 。 

河北 大 海 辽 自然 保护 区 地 处 暖 温带 的 中 山地 区 ,属于 东亚 大 陆 性 季风 气候 区 , 受 地 理 位置 影 响 ,该 区 形成 
的 地 带 性 植被 为 暧 温带 落叶 阔 叶 林 。 山 地 草 甸 分 布 于 2000m 以 上 , 主要 成 分 为 多 年 生 杂 草 ,主要 包括 蔓草 属 
(Carex) .早熟 不 属 ( Poa) 、 路 聚 属 (Papaver) 、 风 毛 策 属 (Saussurea) 金莲 伦 属 (7Trollius) 、 银 莲花 属 (Anemone)、 
马 先 贰 属 ( Pedicularis ) .梅花 草 属 ( Parnassia) 等 


2 数据 来 源 与 研究 方法 


2.1 数据 来 源 
数据 主要 包括 大 海 陀 自 然 保 护 区 近 30 年 的 TM 遥感 影像 数据 ;赤城 县 1980 一 2015 年 的 气候 数据 ,包括 
年 平均 气温 ,7 月 份 平均 气温 、1 月 份 平均 气温 .年 降水 量 .6 一 8 月 平均 降水 量 等 多 个 气候 指标 ;大 海 陀 自 然 保 
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护 区 及 周边 地 区 的 地 形 数据 等 '“ 。 

自然 保护 区 遥感 影像 资料 从 地 理 空 间 数据 云 ( www.gscloud.cn ) 中 下 载 获得 ,使 用 LANDSAT 系列 数据 ,该 
系列 数据 每 16 天 循环 扫描 一 次 “。 具 体 选择 的 遥感 数据 见 表 1。 

地 形 数 据 来 地 理 空间 数据 云 的 CDEMDEM-30M 分 辨 率 数字 高 程 数 据 '*。 

近 36 年 的 温度 和 降水 数据 均 采 用 国家 气象 局 发 布 的 数据 ,温度 和 降水 数据 取 自 1980 一 2015 年 逐年 的 
1km 栅 格 数据 。 


表 1 遥感 数据 基本 信息 


Table 1 Basic information of remote sensing data 


年 份 Year 日 期 (月 -日 )Date(Month-day) 遥感 类 型 Type 数据 标识 Data identification 
1987 9 一 17 Landsat4- 5 TM LT51240321987260BJC00 
1990 9 一 18 Landsat4- 5 TM LT51230321990261HAJ00 
1994 9 一 13 Landsat4- 5 TM LT51230321994256BJC00 
1995 9 一 16 Landsat4- 5 TM LT51230321995259HAJ00 
2001 8 一 31 Landsat4- 5 TM ET51230322001243BJC00 
2004 8—30 Landsat4- 5 T™ LT51240322004243BJCO1 
2006 9 一 21 Landsat4- 5 TM LT51240322006264IKR00 
2007 9—24 Landsat4- 5 T™ LT51240322007267IKR00 
2009 9 一 22 Landsat4- 5 TM LT51230322009265IKRO0 
2010 9 一 25 Léandsat4=5 TM LT51230322010268IKRO0 
2013 9—24 Landsat8 OLI_TRIS LC81240322013267LGN00 
2015 9—14 Landsat8 OL TRIS LC81240322015257LGN00 

2.2 ”数据 人 处理 


2.2.1 归 一 化 差 值 植被 指数 的 计算 

归 一 化 差 值 植被 指数 (Normalized Diffeiepce Vegetation Index, NDVI) 是 反映 植被 禾 新 的 一 个 重要 指 2 导 
时 间 序 列 的 变化 对 应 着 植被 的 生长 和 变化 “| ,NDVI 是 遥感 估算 植被 覆盖 度 研究 中 最 常用 的 植被 指数 ,为 植 
物 生 物 生 长 状态 及 植被 空间 分 布 密度 的 最 佳 指 示 因 子 ,与 植被 分 布 密度 呈 线 性 相关 ,长 期 用 作 监 测 植被 变化 
情况 “ 。NDVI 对 土壤 背景 的 变化 较为 敏感 ,检测 植被 覆盖 度 的 幅度 较 宽 , 当 植 被 覆盖 度 小 于 15% 时 ,土壤 背 
景 与 植被 可 以 区 分 开 ; 当 植 被 覆盖 度 为 23% 一 80% 时 ,NDYVI 值 随 植被 覆盖 度 的 增 大 呈 线 性 增加 ; 当 植 被 覆盖 
度 大 于 80% 时 ,NDVI 的 检测 能 力 逐 步 下 降 ! “7 。 

根据 实地 调查 结果 , 章 铭 主要 位 于 大 海 陀 自然 保护 区 海拔 2000m 以 上 的 地 区 ,通过 GDEMDEM-30M 分 辨 
率 数 字 高 程 数据 提取 该 范围 矢量 文件 ,并 裁剪 计算 平均 NDVI 值 ” 。 

NIR-R 
NIR+R au 

式 中 ,RR 指 红 色 波 段 ,NIR 指 红 外 波段 ,其 中 在 Landsat4-5 TM 分 别 为 Band3 、Band4; 在 Landsat 8 OLI_TRIS 分 
别 为 Band4、Band5。 由 于 季节 对 NDVI 指数 的 影响 ,各 年 度 选 用 相同 的 月 份 用 于 NDVI 指数 计算 ,根据 遥感 影 
像 云 量 . 草 旬 长势 各 月 份 间 长 势 情况 ,选取 8 月 底 至 9 月 下 旬 的 遥感 数据 ,共计 12a. 9 。 

统计 大 海 陀 自 然 保 护 区 海拔 2000m 以 上 地 区 的 平均 NDVI 值 ,采用 Office Excel 与 SPSS 20.0 统计 分 析 软 
件 进 行 数据 分 析 , 然 后 利用 滑动 平均 法 对 数据 进行 处 理 , 得 到 NDVI 值 的 变化 趋势 。NDYVI 取 值 范围 为 -1 一 1， 
其 中 NDVI<0 说 明 无 植被 覆盖 ,NDVI>0 时 有 植被 覆盖 , 且 值 越 大 ,植被 履 盖 度 越 高 ,植被 生长 状况 越 好 。 
2.2.2 气候 变化 数据 整理 与 分 析 

运用 ArcGIS 将 收集 到 的 温度 和 降水 数据 转换 成 面 数据 ,采用 0ffice Excel 与 SPSS 20.0 统计 分 析 软 件 进 
行 单 年 区 域 加 权 平 均 和 多 年 滑动 平均 处 理 , 得 出 温度 和 降水 量 的 变化 趋势 ;运用 基于 秩 的 非 参 数 Mann- 


NDVI= 
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Kendall 统计 检验 方法 (以 下 简称 M-K 检验 法 ) 计 算 温 度 和 降水 量变 化 趋势 是 否 明 显 。 

(1) 清 动 平均 法 

滑动 平均 法 能 够 精确 地 表示 测量 结果 ,抑制 随机 误差 产生 的 影响 。 滑 动 平均 法 即 是 对 非 平稳 的 数据 ,在 
适当 的 小 区 间 上 视 为 接近 平稳 ,从 而 计算 局 部 的 平均 ,以 减 小 随机 起 伏 误 差 , 以 此 类 推 , 沿 全 长 Y 个 数据 逐一 
在 小 区 间 上 进行 不 断 的 局 部 平局 , 即 可 得 出 较 平 滑 的 结果 -” 。 滑 动 平 均 法 常 被 用 于 气候 数据 的 分 析 , 根 据 前 
人 研究 经 验 “” ,本 文选 取 5 年 滑动 平均 法 “ 。 

请 动 平 均 法 为 沿 全 长 N 个 数据, 不断 逐一 滑动 取 m 个 相 邻 数据 进行 加 权 平 均 运 算 ,得 出 平滑 数据 ,其 一 
般 算 式 为 : 


p 
Pe lr 2 (2 ) 
7 二 4 


式 中 ,w, 为 权 系数 ， 立 w = 1; 且 79 为 小 于 疾 的 任 一 正 整数 , 且 Pro+1=mla 

运用 统计 学 方法 一 元 线性 方程 定量 描述 气候 因子 的 回归 变化 趋势, 即 以 年 代为 栖 间 因子 ;计算 保护 区 的 
年 均 温 .1 月 平均 气温 .7 月 平均 气温 、 年 降水 量 .6 一 8 月 平均 降水 量 的 滑动 平均 数据 的 线性 趋势 。 以 年 份 为 
自 变量 ,要 素 (气温 、 降 水量 ) 分 别 为 因 变量 ,建立 一 元 回归 方程 对 要 素 和 年 份 数据 进行 线性 回归 分 析 , 计 算 直 
线 的 拟 合 斜率 5 和 相关 系数 尺 值 。 用 蕊 为 第 ;年 要 素 (气温 、 降 水 量 闪 用 爷 表示 六 所 对 应 的 时 间 , 建 立 一 元 线 


Xi=atbt(i=1,2,5n) (3) 
| 一 SU AAS ;) 
本 i=1 让 1 i=1 
7 二 人 ( Dy2 
#1 B=1 


Yt ZY) (Y; se 
R= z -外 (35) 


> 


式 中 ,oa 为 截 距 ,为 斜率 (变化 趋势 ) ,a 和 bb 运用 最 小 二 乘法 估计 。 当 和 斜率 45>0 时 ,说 明 该 要 素 变化 随 年 度 呈 
上 升 趋势 ; 若 b<0 时 ,说 明 该 要 素 变化 随 年 度 呈 下 降 趋势 “ 。 

(2) M-K 检验 法 

该 方法 常用 于 水 质 ,人 径流 量 、 温 度 、 降 水 等 时 间 序 列 变 化 趋势 的 显著 性 检验 ,主要 是 通过 计算 统计 量 7, 方 
差 o, 和 标准 化 变量 民 ; 来 判断 序列 趋势 是 否 显著 : 


b (4) 


M = 一 (0) 
OO, 
45 
| C7 
MA 
9NOV- DT ， 


式 中 ,5 为 序列 中 数值 间 i<i 时 ,< 出 现 的 次 数 ;其 中 Jj 为 年 份 ,XX 对 分 别 为 i 年 的 数据 ;NN 为 序列 长 度 ; 
当 取 w= 0.05 的 显著 水 平时 , 耕 1M1>1.96 ,说 明 变 化 显著 ; 当 取 w=0.01 的 显著 水 平时 , 奋 1M1>2.58 ,说 明 变 化 
极 显著 ;W 值 为 正 ,表明 具有 上 升 或 增加 趋势 ,W 值 为 负 , 则 意味 着 下 降 或 减少 的 趋势 。 
2.2.3 影响 草 甸 变化 的 因素 分 析 

地 理 要 素 之 间 做 相关 分 析 主 要 目的 是 揭示 各 要 素 之 间 的 密切 程度 ,而 这 种 密切 程度 的 测定 , 则 要 通过 计 
算 它 们 之 间 的 相关 系数 并 检验 来 进行 。 
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(1) 相 关 性 分 析 

相关 分 析 是 考量 两 个 连续 变量 之 间 的 相关 关系 ,实际 对 任何 类 型 的 变量 ,都 可 以 运用 相关 的 指标 分 析 两 
者 之 间 的 相关 关系 , 它 能 以 数字 方式 准确 地 描述 变量 间 的 线性 相关 程度 。 

Spearman 秩 相关 系数 是 基于 各 因子 间 秩 的 大 小 ,属于 非 参数 检验 ,其 表达 式 为 : 

6 
(i) =1- -全 全 (9) 

式 中 ,r (7 为 Spearman 秩 相 关系 数 ;n 为 因子 总 年 数 ;d, = (Ce ,Xj 和 寻 分 别 为 因子 和 因子 7 在 左 年 份 
中 的 重要 值 ” 。 

(2) 偏 相关 分 析 

偏 相 关系 数 是 在 对 其 他 变量 影响 进行 控制 的 条 件 下 ,衡量 多 个 变量 中 某 两 个 变量 之 间 的 线性 相关 程度 的 
指标 。 通 过 偏 相 关系 数 与 相关 系数 的 比较 ,来 确定 两 个 变量 之 间 的 内 在 线性 联系 会 更 真实 更 可 靠 ”。 偏 相 
关 分 析 的 计算 方法 是 : 先 计 算 3 个 因子 间 的 相关 系数 ,然后 通过 3 个 人 简单 相关 系数 来 计算 仿 相 关系 数 。 其 公 
式 为 : 


Tij(k) 一 7 
ya Ee V1 加 0 
式 中 ,r,(h) 就 是 在 控制 了 因子 影响 所 计算 的 因子 i 和 因子 j 之 间 的 偏 相关 系 烤 yr .rj .rj 分 别 是 i .43 个 因 
子 之 间 的 两 两 简单 相关 系数 。 


(10) 


3 结果 分 析 


3.1 草 甸 植被 NDVI 变化 07 
讨论 NDVI 变化 特征 时 ,生长 值 的 峰值 作为 备用 的 ” 0 
特征 量 之 一 ,能 很 好 地 反映 植被 生长 发 育 程 度 和 未 来 发 _ 


展 趋势 等 , 故 选取 8 月 底 至 9 月 下 凶 秆 宜 此 长 天 沉 ” 急 。 
较 好 时 有 段 的 遥感 影像 进行 分 析 '“( 表 了 )， 0.2 
由 大 海 陀 自然 保护 区 1987 一 2015 年 间 选 取 的 12 0.1 
年 的 草 甸 植被 NDVI 指数 变化 趋势 图 (图 /1 ) 可 以 看 出 ， "人 SASa 和 Sav 
草包 植被 NDVI 在 1987 一 2004 年 呈现 小 幅 持 续 增 加 的 四 
趋势 ,在 2004 一 2015 年 间 旦 下 降 趋势 ,NDVI 值 在 2004 
年 达到 最 高 值 058- 随后 逐渐 下 降 ,在 2013 年 达到 了 最 图 1 大 海 陀 自然 保护 区 草 多 NDVI 指 关 
低 值 0 23 [0] 6 Fig.1 The NDVI index of Dahaituo Nature Reserve meadow 


3.2 ”区 域 气候 变化 
3.2.1 降水 量变 化 

为 了 获取 研究 期 (1980 一 2015 年 ) 内 大 海 陀 自然 保护 区 的 降水 量变 化 ,采用 统计 学 方法 计算 了 该 时 期 的 
年 降水 量 及 6 一 8 月 平均 降水 量 ,结果 如 图 2、 图 3 所 示 '。 

从 图 2 .图 3 中 可 以 看 出 ,大 海 陀 自然 保护 区 草 馈 区 域 近 36 年 的 年 降水 量变 化 整体 表现 为 波动 循环 的 趋 
势 , 旦 总 体 略 有 上 升 , 其间 伴随 有 小 幅 波 动 。 其 中 年 降水 量 2007 年 最 低 ,为 619.48mm; 2000 年 最 高 , 达 
1356.35mm;6 一 8 月 平均 降水 量 2007 年 最 低 ,为 102.46mm;1997 年 最 高 , 达 230.79mm ” 。 

从 年 代 变化 来 看 , 1980 一 1991 年 降水 量 平均 为 915.04mm,1992 一 2003 年 降水 量 平均 为 950.92mm， 
2004 一 2015 年 降水 量 平均 为 912.42mm ,长 期 变化 的 趋势 并 不 明显 !9 。 

经 M-K 检验 法 检验 ,M( 年 降水 量 )= 0.74<1.94,M(6 一 8 月 平均 降水 量 )= 1.20<1.94, 均 未 达到 显著 水 
平 ,也 说 明了 大 海 陀 自 然 保 护 区 年 降水 量 和 6 一 8 月 平均 降水 量 在 近 36 年 的 时 间 序 列 上 变化 不 显著 “| 。 
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图 2 1980 一 2015 年 大 海 陀 自然 保护 区 6 一 8 月 平均 降水 量变 化 趋势 


Fig.2 The change trend of average precipitation in June to August in Dahaituo Nature Reserve meadow during 1980 to 2015 
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图 3”1980 一 201S 年 大 海 陀 自 然 保护 区 年 降水 量变 化 趋势 


Fig.3 The change trend ofannual precipitation in Dahaituo Nature Reserve meadow during 1980 to 2015 


3.2.2 人 气温 变化 

(1) 年 均 温 变化 趋势 

1980 一 2015 年 大 海 陛 自然 保护 区 草 多 气温 变化 趋势 如 图 4 所 示 。 近 36 年 该 地 区 的 平均 气温 为 6.47C， 
年 平均 增 温 速 率 为 0.55%/A10a。 从 年 代 变 化 来 看 ,1980 一 1991 年 均 温 为 5.92%C ,1992 一 2003 年 均 温 为 6.14%C ， 
2004 一 2015 年 均 温 为 7.34%C ;其 中 虽然 在 1980 一 1991 年 间 有 阶段 性 下 降 的 趋势 ,但 从 总 体 来 看 仍 呈 阶段 性 上 
升 的 趋势 ;期 间 的 最 高 年 均 温 为 7.95% ,出 现在 2007 年 ;最 低 年 均 温 为 5.02%C ,出 现在 1980 年 '*1。 

经 M-K 检 验 法 检验 ,MM( 年 均 温 )= 4.54>2.58 ,到 达 了 极 显 著 水 平 (P<0.01) ,说 明了 该 地 区 的 气温 上 升 趋 
势 极 显著 ,与 全 球 变 暖 的 趋势 较为 一 致 "” 。 

k9) 1 月 7 月 平均 气温 变化 趋势 

图 5、 图 6 为 大 海 陀 自然 保护 区 草 旬 1 工 月 7 月 的 平均 气温 变化 趋势 。 该 地 区 1 月 .7 月 的 月 平均 气温 均 呈 
现 出 了 上 升 的 趋势 ,1 月 平均 气温 的 上 升 趋势 明显 高 于 7 月 平均 气温 的 上 升 趋势 ,1 月 的 气温 增 温 速 率 为 
0.76% /10a;7 月 气温 增 温 速 率 为 0.27%C7Z10a 9 。 

经 M-K 检验 法 检验 ,M(1 月 均 温 )= 2.10>1.94,M(7 月 均 温 )= 0.65<1.94, 其 中 1 月 平均 气温 的 变化 到 达 
了 显著 水 平 (P<0.05) ,说 明了 该 地 区 的 气温 上 升 趋势 主要 由 以 1 月 为 主 的 冬季 升 高 而 引起 。 


http://www.ecologica.cn 


7 期 宋 超 


等 : 近 30 年 河北 大 海 陀 自然 保护 区 山地 草 多 植被 (NDVI) 变化 及 其 对 气候 的 响应 


2553 


年 均 气 温 
Annual average temperature/°C 


1980 1983 1986 1989 


1992 


1995 


一 人 一 年 均 温 
一 各 一 5a 滑动 平均 
一 一 线性 变化 趋势 


1998 2010 2013 2015 


年 份 Year 


2001 2004 2007 


图 4 1980 一 2015 年 大 海 陀 自然 保护 区 年 均 温 变 化 趋势 


Fig.4 The change trend of average annual temperature in Dahaituo Nature Reserve meadow during 1980 to 2015 
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图 5 1980 一 2015 年 大 海 陀 自然 保护 区 1 月 份 平 均 气 温 变 化 趋势 


Fig.5 The change trend of average temperature in January in Dahaituo Nature Reserve meadow during 1980 to 2015 


3.3 ”草包 植被 NDVT 与 气候 变化 的 关系 


3.3.1 草 甸 植被 NDVI 与 气候 变化 的 相关 分 析 
分 别 将 草 甸 植被 NDVI 值 与 降水 量 温度 等 环境 因子 作 相 关 分 析 ,结果 见 表 2。 
表 2 草 旬 植被 NDVI 与 环境 因子 的 相关 性 分 析 


Table 2 The correlation analysis between NDVI of the meadow and environmental factors 


环境 因子 Spearman 秩 相 关系 数 
Eiivironmental factor Spearman's rank correlation coefficient 

年 均 温 Average annual temperature -0.664 0.018* 
7 月 均 温 Average temperature in July -0.168 0.602 

1 月 均 温 Average temperature in January -0.769 0.003 
年 降水 量 Annual precipitation 0.084 0.795 

6 一 8 月 平均 降水 量 Average precipitation (June to August) 0.098 0.762 


* P<0.05; * P<0.01 


http://www.ecologica.cn 


2554 生 态 学 报 38 卷 


7 月 平均 温 


The average temperature in July/°C 
BS 
| 
> 
省 
| 
{ 
有 
Re 
| 


19 
1980 1983 1986 1989 1992 1995 1998 2001 2004 2007 2010 2013 2015 


年 份 Year 


图 6 1980 一 2015 年 大 海 陀 自然 保护 区 7 月 份 平均 气温 变化 趋势 
Fig.6 The change trend of average temperature in July in Dahaituo Nature Reserve meadow during 1980 to 2015 


在 Spearman 秩 相 关 分 析 中 可 以 看 出 , 草 甸 植 被 NDVI 与 温度 均 晨 人 负 相 关 , 且 与 年 均 温 及 1 月 均 温 达到 了 
显著 相关 水 平 (P<0.05); 且 相 比 之 下 ,1 月 均 温 的 显著 性 指数 达到 了 =0.769。 说 明 当 年 的 年 均 温 越 低 , 草 多 植 
被 NDVI 指数 越 高 , 草 甸 植 被 生长 状况 越 好 , 且 主 要 受 1 月 均 漫 的 影响 , 受 7 月 均 温 的 影响 较 小 1 。 

草 多 植被 NDVI 与 年 降水 量 .6 一 8 月 平均 降水 量 均 呈正 相关 ,但 均 不 显著 
3.3.2 ” 草 旬 植被 NDVI 与 气候 变化 的 偏 相关 分 析 

分 别 将 代表 草 馈 植被 NDVI 与 降水 量 、 年 均 温 等 环境 因子 的 变化 作 偏 相关 分 析 , 分 别 控制 变量 年 均 温 和 
年 降水 量 , 对 草 甸 植被 NDVI 变化 进行 分 析 ( 表 3), 可 以 看 出 , 草 甸 植被 NDVI 与 年 均 温 的 偏 相 关系 数 为 
-0.738 ,通过 P<0.01 的 显著 性 水 平 检验 ,上 且 负 相关 程度 大 于 未 作 偏 相关 分 析 的 相关 系数 ,说 明 排除 了 年 降水 
量 影响 后 ,年 均 温 与 草 馈 植被 NDVI 的 相关 程度 更 加 明显 ; 草 多 植被 NDVI 与 年 降水 量变 化 的 偏 相 关系 数 为 
-0.249 ,未 通过 显著 性 水 平 检验 ,说明 草 甸 植被 NDVI 与 同期 降水 量 属于 不 显著 的 负 相 关 '| 。 


表 3 章 匈 植被 NDVI 与 环境 因子 偏 相关 分 析 


Table 3 The partial correlation analysis between NDVI of the meadow and environmental factors 


控制 变量 Control variable 偏 相 关 > Partial correlation P 显著 性 Significant 
年 均 温 Average annual temperature —0.249 0.461 不 显著 
年 降水 量 Annual precipitation 一 0.738 0.009 显著 


4 结论 与 讨论 


4:1 结论 

(1) 大 海 陀 自然 保护 区 草 甸 植被 NDVI 值 呈 先 上 升 后 下 降 趋 势 ,在 2004 年 达到 最 高 值 0.58 ,在 2013 年 降 
到 最 低 值 0.23。 

(2) 大 海 陀 自然 保护 区 草 旬 区 域 的 降水 量 .6 一 8 月 平均 降水 量 呈 波动 循 坏 趋势 ,总 体 略 有 上 升 ,但 均 未 达 
到 显著 水 平 。 该 地 区 年 均 降 水 量 为 915.04mm ,在 2000 年 年 降水 量 达 到 最 高 值 1356.35mm ,在 2007 年 年 降水 
量 降 到 最 低 值 201.46mm 。 

通过 Spearman 秩 相 关 分 析 , 草 旬 植 被 NDVI 与 年 降水 量 .6 一 8 月 平均 降水 量 呈 正 相 关 ,但 均 不 显著 。 通 
过 偏 相关 分 析 , 草 旬 植 被 NDVI 与 年 降水 量变 化 的 偏 相 关系 数 为 -0.249 ,未 通过 显 若 性 检验 ,为 不 显著 负 相 关 
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(3) 大 海 陀 自然 保护 区 草 甸 区域 的 年 均 温 1 月 平均 气温 7 月 平均 气温 均 呈 上 升 趋势 ,其 中 年 均 温 达 到 
极 显著 水 平 ,1 月 平均 气温 达到 显著 水 平 。 该 地 区 年 均 温 为 6.47%C ,在 2007 年 ,年 均 温 达到 最 高 值 7.95%C ,在 
1980 年 年 均 温 达到 最 低 值 5.02%C 。 

通过 Spearman 秩 相 关 分 析 , 草 甸 植 被 NDVI 与 温度 均 呈 负 相 关 , 且 与 年 均 温 及 1 月 平均 气温 为 显著 负 相 
关 。 说 明 当 年 的 年 均 温 越 低 , 草 旬 植 被 NDVI 指数 越 高 , 草 旬 植 被 生长 状况 越 好 , 且 主 要 受 1 月 均 温 的 影响 ， 
受 7 月 均 温 的 影响 较 小 。 通 过 偏 相关 分 析 , 草 甸 植 被 NDVI 与 年 均 温 的 偏 相关 系数 为 -0.738 ,通过 显著 性 检 
验 ,说明 排除 了 年 降水 量 影响 后 ,年 均 温 与 草 匈 植被 NDVI 的 相关 程度 更 加 明显 。 

4.2 ”讨论 

在 山地 自然 生态 系统 中 ,地 形 因 子 能 影响 温度 和 降水 的 再 分 配 ,一 般 而 言 , 随 着 海拔 高 度 的 上 升 ;气温 还 
渐 降 低 ,降水 量 和 相对 湿度 呈 增 加 趋势 。 植 被 生长 在 低 海拔 地 区 主要 受到 水 分 条 件 的 限制 ,降水 量 的 增加 将 
促进 植被 的 生长 ,在 高 海拔 区 域 , 植 被 生长 会 受到 热量 条 件 的 限制 ,温度 上 升 ,或 者 温度 上 升 和 降水 量 增加 共 
同 改善 了 植被 的 生长 条 件 , 有 利于 植被 生长 ”*” 。 根据 以 上 研究 中 对 草 甸 所 在 区 域 海拔 的 划分 ,本 人 研究 中 大 
海 陀 自然 保护 区 草 甸 分布 区 属于 低 海拔 地 区 , 近 30 年 来 ,降水 没有 显著 变化 , 草 馈 植被 NDVI 与 年 降水 量 也 
没有 显著 的 相关 性 ,但 却 与 年 均 温 表现 为 负 相 关 性 ,在 降水 量 不 变 的 情况 下 , 仪 温度 升 高 s 反 而 对 草 甸 植被 的 
生长 具有 抑制 作用 , 即 年 均 温 升 高 ,NDVI 指数 降低 。 有 学 者 在 对 我 国 东 北 多 年 冻 土 区 植被 NDVI 的 研究 后 认 
为 ,伴随 着 气温 的 显著 升 高 和 降水 量 下 降 , 植 被 NDVI 显著 下 降 ,降水 量 是 影响 植被 NDVI 的 主要 因子 ,这 
与 本 研究 得 出 的 结果 类 似 “ 。 

姚 玉 壁 等 人 在 对 甘肃 玛 曲 县 亚 高 山 草 旬 人 研究 后 ;发 现 山地 草 和 甸 类 草地 植被 从 返青 到 分 药 期 ,气温 对 草地 
牧草 产量 形成 正 效 应 ,拔节 到 抽 惩 期 牧草 产量 对 温度 变化 敏感 ;开花 到 籽 实 成 熟 期 ,热量 对 牧草 产量 的 影响 由 
正 效应 转向 负 效 应 ,在 籽 实 成 熟 期 牧 革 产量 形成 对 气温 变化 进入 第 二 个 敏感 期 ; 除 籽 实 成 熟 期 降水 量 对 牧草 
产量 形成 为 负 相 关外 ,其余 时 段 降水 量 对 牧草 产量 形成 均 为 正 效 应 "” ,本 研究 选取 8 月 底 至 9 月 下 旬 草 甸 遥 
感 影像 资料 ,根据 姚 玉 壁 的 草 旬 植 被 生长 期 分 类 ;大 海 陀 自然 保护 区 草 甸 在 8 月 底 至 9 月 下 旬 为 籽 实 期 至 黄 
村 期 ,温度 对 牧草 产量 影响 为 负 效 应 ,与 本 研究 得 出 的 结论 一 致 “ 。 

在 对 我 国 东 部 草 甸 所 在 地 区 的 植被 覆盖 变化 研究 中 , 张 建 亮 等 人 对 长 白山 地 区 的 的 研究 中 采用 了 7 一 8 
月 份 NDVI1, 发 现 植被 变化 受气 温 的 影响 程度 均 强 于 降雨 ,温度 的 波动 对 植被 的 影响 更 为 敏感 。 黄 晓 霞 
等 人 对 小 五 台 亚 高 山 草 多 的 研究 中 ?从 景观 格局 的 角度 研究 发 现 ,制约 草 多 群落 分 布 的 最 为 重要 的 条 件 是 热 
量 状况 ,其 次 是 水 分 条 件 的 影响 。 这 些 与 本 研究 得 出 的 结论 一 致 。 

大 海 陀 自然 保护 区 草 多 的 植被 变化 除了 受降 雨量 和 温度 等 气象 因子 影响 外 ,与 旅游 活动 等 人 为 因素 有 一 
定 关 系 。 随 着 当地 政府 对 保护 区 的 保护 力度 加 大 , 草 甸 植 被 逐步 改善 ,大 海 陀 自 然 保护 区 于 2003 年 被 评 为 
家 级 自然 保护 区 同时, 当地 政府 加 大 了 保护 区 的 投资 开发 力度 ,大 海 陀 逐渐 成 为 北京 周边 旅游 的 热门 地 区 
之 一 , 草 甸 成 为 游人 活动 的 集中 区 域  。 但 由 于 地 处 相对 偏僻 ,社会 经 济 发 展 比 较 落 后 ,对 旅游 的 管理 落后 ， 
缺乏 科学 规划 和 专业 人 员 ,可 能 成 为 导致 草 甸 植 被 生长 状况 出 现 逐 步 下 降 的 因素 之 一 。 
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